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We examine the structural and dynamic properties of confined binary hard-sphere mixtures designed to 
mimic realizable colloidal thin films. Using computer simulations, governed by either Newtonian or over¬ 
damped Langevin dynamics, together with other techniques including a Fokker-Planck equation-based 
method, we measure the position-dependent and average diffusivities of particles along structurally isotropic 
and inhomogeneous dimensions of the fluids. At moderate packing fractions, local single-particle diffusivities 
normal to the direction of confinement are higher in regions of high total packing fraction; however, these 
trends are reversed as the film is supercooled at denser average packings. Auxiliary short-time measure¬ 
ments of particle displacements mirror data obtained for experimental supercooled colloidal systems. We 
find that average dynamics can be approximately predicted based on the distribution of available space for 
particle insertion across orders of magnitude in diffusivity regardless of the governing microscopic dynamics. 

PACS numbers: Valid PACS appear here 


Confined fluids exhibit inhomogeneous structural and 
relaxation properties, which are general features of ma¬ 
terials subjected to position-dependent external fields. 
Because confined fluids emerge in a diverse array of nat¬ 
ural and technological contexts (e.g., water in biological 
media, polymer thin films, etc.), considerable attention 
has been directed at understanding how their static and 
dynamic properties relate to bulk fluid physics observed 
under similar conditions. As a result, the static prop¬ 
erties of confined fluids, such as local one-body den¬ 
sity p{z), are now well-understood in terms of physi¬ 
cal intuition (e.g., emergence of particle layering near 
boundaries to relieve packing frustration^’^) and can 
be predicted using microscopic approaches like density 
functional theory^’^. However, much less is understood 
about what controls the dynamics of inhomogeneous flu¬ 
ids, and only recently have efforts broadened to include 
developing theories^^^^ and other tools^^“^® for charac¬ 
terizing particle dynamics both on a spatially-averaged 
basis and as a function of position. 

Given the difficulty of applying first principles to un¬ 
derstand the dynamics of such systems, progress has 
been made by virtue of use pragmatic approaches, e.g., 
application and testing of semiempirical, quasi-universal 
scaling laws that relate transport coefficients of inter¬ 
est to static properties^^^^"^. To wit, it has been shown 
that single-particle diffusivities, relaxation times, and 
viscosities along structurally-invariant (i.e., isotropic) 
dimensions of simple confined fluids can be predicted 
based on knowledge of how dynamic properties of the 
bulk fluid relate to static quantities including excess en¬ 
tropy (relative to the ideal gas) and fractional avail¬ 
able space exp{c^^^} (or insertion probability po), which 
characterize short-range static correlations and particle 
packings, respectively^®^^^. 

In this spirit, one might expect that local particle mo¬ 


bility in an inhomogeneous fluid should similarly corre¬ 
late with position-dependent static properties; in other 
words, the way particles navigate through the inhomo¬ 
geneous environment might be encoded in the physics of 
motion observed in a bulk, homogeneous fluid. However, 
the validity of such a connection has yet to be carefully 
and systematically evaluated. Despite providing other 
important insights, previous investigations directly mea¬ 
suring inhomogeneous dynamics have studied a variety 
of fluids governed by disparate interactions, external 
fields, and conditions, and they have also used differ¬ 
ent protocols to characterize the dynamics®’^^’^^’^^’^®“‘^°. 
As a result, even fundamental questions related to con¬ 
fined (and more generally inhomogeneous) fluids remain 
open: Do local and average correlations between parti¬ 
cle mobility and structure universally reflect bulk be¬ 
haviors? Do new structure-mobility relations emerge as 
inhomogeneous fluids are supercooled toward glass tran¬ 
sitions? And does the choice of microscopic dynamics 
affect these qualitative trends? 

As a step toward addressing these questions, we ex¬ 
amine computer simulations of bulk and confined binary 
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FIG. 1. (color online). Illustration of confined binary mix¬ 
ture comprising small (blue) and large (red) particles. 
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mixtures of small (sm) and large (Ig) hard spheres (HS) 
approximated by a steeply-repulsive Weeks-Chandler- 
Andersen (WCA) pair potential^^ between particles i 
and j, adapted for multiple particle diameters: '•Pi,j(r) = 
4e([crsm/(r + A)]^® - [CTsm/(?’ + A)]^'^) + £ for r < 

- A) and (fij{r) = 0 for r > (2^/24^^^ _ 
where e is the characteristic energy scale; r is the in¬ 
terparticle separation; a denotes particle diameter; and 
A = cTsjjj — (l/2){(Ji + aj). The binary mixtures are com¬ 
posed of spheres with size ratio CTig/cTsm = 1-3, volume- 
proportional masses mig/msm = (cig/CTsm)^) and com¬ 
position defined by the fraction of small particles Xsm = 
0.75. These parameters mimic colloidal mixtures inves¬ 
tigated in recent experiments^^’®®. Below, we implicitly 
non-dimensionalize quantities via appropriate combina¬ 
tions of the characteristic lengthscale CTsm and energy 
scale e = fcaT, where is Boltzmann’s constant, and 
T is temperature. 

Periodic boundary conditions are applied in all direc¬ 
tions for the bulk systems, while for the confined sys¬ 
tems (see Fig. 1), particles are situated in slit-pores 
of size H = 5 between two reflective walls placed at 
z = ±i7/2, with periodic boundary conditions applied 
in the x- and y-directions. The wall-particle interac¬ 
tions are analogous to the hard-sphere-like interactions 
between particles, but defined such that the center of 
particle i can access —{H — cri)/2 ^ z < {H — CTi)/2. 
Spatially averaged packing fractions are given by = 
(7r/6)/0*[xsm + 0 ’fg(l - Xsm)], where p* = {Nsm + Nig)/V 
is the combined number density of both species and V 
is volume. Here, the values for the confined fluids 
are dehned by the total (surface- rather than center- 
accessible) slit pore volume. We generate particle tra¬ 
jectories governed by either conventional molecular dy¬ 
namics (MD) or Brownian dynamics (BD) (i.e., over¬ 
damped Langevin ignoring hydrodynamic interactions) 
using GROMACS 4.5.5^^ with implementation details 
are provided in the Supplemental Material (SM)'^®. 

To characterize particle motions, we calculate mean- 
squared displacements (MSDs) and diffusivities in the 
structurally isotropic and inhomogeneous directions of 
the bulk and confined systems. Average diffusivities 
£)avg bulk systems and parallel to the walls in 

the confined systems characterize motions in isotropic 
directions, and are derived by fitting the long-time be¬ 
havior of the MSD of all the particles to the Einstein 
relation (Ar^) = 2dDAt. In the bulk (confined) case, 
(Ar^) is the MSD in the x-, y-, and z-directions {x- and 
y-directions) over lag-time At and dimensionality d = 3 
(d=2). 

Diffusivities in the inhomogeneous z-direction of the 
confined pores are position-dependent and cannot be 
calculated via the Einstein relation because particles are 
subjected to locally non-cancelling potentials of mean 
force"^^. Particle displacements along the z-coordinate 
are instead accurately described®® by the ID Fokker- 
Planck (FP) equation 

^ = ( 1 ) 



_^^^^_ r , _3 
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FIG. 2. (color online). Local total packing fractions (f>{z) 
(top) and local diffusivities in the z-direction D^{z) of small 
particles calculated from MD simulations (middle) and BD 
simulations (bottom) of pore size H = 5 and average total 
packing fractions = 0.20, 0.30, 0.35, 0.40, 0.45, 0.48, 
0.50, 0.51, and 0.52. The BD profiles are normalized by 
the infinite dilution diffusivity D^iyz) profiles for the 

large particles exhibit shapes in line with the small-particle 
profiles, as shown in the SM®®. 


where Dz(z) are position-dependent diffusivities. Here, 
G{z,to + At\z',to) is the Markovian propagator char¬ 
acterizing temporal single-particle displacements given 
the potential of mean force F{z) = — ln{p(z)} -|- C, 
where C is an arbitrary constant. To obtain D^{z) 
from simulation data, we use a mean-first passage times 
(MFPT) method®®’®®’®® applied to the steady-state (i.e., 
dGjdt = 0) limit of the FP equation, which is known to 
provide equivalent information compared to alternative 
FP treatments®®’®®. Additional implementation details 
are discussed elsewhere®®’®®. 

We begin our discussion by considering Fig. 2, where 
we compare local total packing fraction 4'{z) and local 
particle diffusivities in the z-direction Dz{z) for con¬ 
fined systems over a wide range of and governed 
by either Newtonian or Brownian microscopic dynamics. 
Here, we use (j){z) because it more economically quan¬ 
tifies the local aggregate packing frustration compared 
to component density profiles p(z), where the latter are 
provided in the SM®®. Remarkably, we find that while 
the packing structure in the confined pores undergoes 
an apparent shift from four to five dense particle layers 
upon increasing the shapes of the Dz{z) profiles 

are qualitatively insensitive to this considerable strue- 
tural rearrangement. Thus, for < 0.40, particles 
diffuse more quickly though densely-packed regions (ex¬ 
cept very close to the walls, where particles slow down 
due to impenetrability), but for > 0.45, particles 
instead move more slowly through densely-packed re- 
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FIG. 3. (color online). Mean-squared displacements (MSD) 
per particle in the x- and y-directions versus lag-times At 
for small (solid lines) and large (dashed lines) particles from 
MD (main) and BD (inset) simulations of pore size H = 5 
and average total packing fractions = 0.20, 0.30, 0.35, 
0.40, 0.45, 0.48, 0.50, 0.51, and 0.52. In (b), only small- 
particle curves are shown for clarity and lag-times have been 
normalized by Symbols in (a) denote times cor¬ 

responding to the prohles in Fig. 4. 


gions. 

This gradual reversal from positive to negative lo¬ 
cal correlations between packing fraction and mobility 
bridges observations based on previous measurements 
of particle dynamics in confined pores-measurements 
that seemingly pointed to inconsistent local trends, but 
where comparisons were also complicated by different 
protocols and dynamic regimes. Mittal et al.^^ mea¬ 
sured local FP-based diffusivities in Newtonian HS sim¬ 
ulations at equilibrium conditions < 0.40) and ob¬ 
served positive correlations between local density p{z) 
(or (^(z)) and D^{z). In contrast, Nugent and co- 
workers^'^’^® experimentally measured short-time MSDs 
along the 2 -coordinate as a function of position for su¬ 
percooled thin films of pseudo-HS colloids, results which 
pointed to negative correlations between local density 
and mobility. While the latter results more intuitively 
correlate to expectations based on bulk HS density 
trends, Mittal et. al. provide a plausible physical basis 
for the observed positive correlations. Specifically, they 
correctly note that higher-density regions in such inho¬ 
mogeneous HS systems also exhibit the greater fraction 
of locally available space for inserting additional parti¬ 
cle centers, i.e., more locally free volume, which might 
correlate with dynamics^’^^’"'®. 

The results in Fig. 2 clearly demonstrate within a 
single framework for measuring dynamics that either 
positive or negative correlations between density and 
diffusivity can be observed in these systems depending 
on whether equilibrium or supercooled conditions are 
being studied. By considering Fig. 2 in conjunction 



FIG. 4. (color online). Local total packing fractions (f>{z) 
and individnal component densities p{z) (top) and MSDs per 
particle in the x- and j/-directions (middle) and 2 -direction 
(bottom) for various lag-times At plotted as a function of 
particle position zo at At = 0. Apart from 4>{z) profiles 
(line-symbols), results for small and large particles are plot¬ 
ted with darker and lighter curves, respectively. Results are 
calculated from MD simulations for pore size iL = 5, where 
left panels show results for = 0.35 and At = 0.05, 0.1, 
0.2, 0.5, 2.0, and 10.0 and right panels for = 0.52 and 
At = 0.1, 1.0, 50.0, and 500.0. Lag-times are also plotted in 
Fig. 3 as symbols. 

with Fig. 3, where the latter shows particles MSDs par¬ 
allel to the confining walls as a function of lag-time At 
and we observe that the reversal in local structure- 
mobility correlations approximately coincides with the 
emergence of plateaus in the MSDs at > 0.45, a 
signature of sub-diffusive “particle caging” characteris¬ 
tic of supercooling^'^ds 

The above results imply that local packing structure 
as measured by (j){z) does not generally correlate in a 
nontrivial way to position-dependent diffusive mobility 
(this is also true of more “microscopic” local static quan¬ 
tities like po{z), as shown in the SM'^®). In turn, given 
that bulk HS fluids exhibit simple negative correlations 
between packing fraction and mobility, it is apparent 
that local static-dynamic correlations in confined-fluid 
systems cannot be naively extrapolated (or predicted) 
from the bulk physics, in agreement with findings for 
more idealized density-varying HS systems®®. Interest¬ 
ingly, for the systems examined here, the choice of mi¬ 
croscopic dynamics had no qualitative impact on the 
shapes of the D^{z) profiles, though recent results"^® in¬ 
dicate this is not generally true of inhomogeneous fluids. 

We next provide data reinforcing the idea that the 
opposing correlations in Fig. 2 between packing frac¬ 
tion and diffusivity at equilibrium versus supercooled 
conditions may also emerge in real colloidal thin films: 
in Fig. 4, we compare position- and time-dependent 
particles MSDs of the confined films against (j){z) and 













































4 


component p(z), results that mirror experimental mea¬ 
surements by Nugent and co-workers (see, e.g., Figs. 
7-8 from^®) for supercooled thin films of pseudo-HS. In 
particular, we show results calculated from MD simu¬ 
lations for = 0.35 and 0.52, which correspond to 
equilibrium and supercooled conditions, respectively. 

Mirroring the experimental findings, for the super¬ 
cooled conditions in Fig. 4, MSDs in the xy-plane are 
insensitive with respect to originating position zq (i.e., 
position at lag-time At = 0) for all At, while MSDs 
in the inhomogeneous ^-direction are negatively corre¬ 
lated with respect to ^(z) (p(z)) for sufficiently short 
At. At longer At, the MSD dependence on zq disap¬ 
pears as particles are no longer generally situated near 
zo- At = 0.35, MSDs in the xy-plane also do 
not vary with zq, but MSDs in the z-direction are in¬ 
stead positively correlated with respect to (j){z) (p(z)) 
for At < 0.2, though the correlations appear to reverse 
at longer At > 0.50 before washing out at long lag- 
times. To our knowledge, no analogous data for real 
non-supercooled thin films for has yet been published. 

Surprisingly, the MSDs in the z-direction at many 
At approximately reflect the Dj,(z) profiles in Fig. 2 
even though the quantified motions are by necessity sub- 
diffusive and accrued when a particle is no longer at zq. 
In turn, only profiles at the very shortest times provide 
information about motions precisely at zq, but these 
are also furthest from the diffusive regime. Nonethe¬ 
less, given that the results for the supercooled system 
in Fig. 3 are consistent with the available experimental 
data, it is plausible that the FP-derived results for diffu¬ 
sive motions-including the positive correlations between 
(j){z) and Dz(z) and their reversal at high be 

observed in real confined colloids. 

Given that the previous results undermine any notion 
of a universal connection between local structure and 
local mobility, it is natural to wonder whether average 
diffusivities in directions parallel and perpendicular to 
the confining walls reflect bulk fluid physics and can be 
predicted based on average static properties. In Fig. 5, 
we address this by comparing average diffusivities DJy® 

and D^"s(z) = D,(z)p(z)dz j ff^%(z)dz for the 
small and large particles from the confined pores against 
curves for bulk mixtures. Here, we plot these dynamic 
quantities against component-specific average insertion 
probabilities pg'"® (or available volumes for insertion), 
which have been shown to provide the most quantita¬ 
tively robust connection between bulk diffusivity and 
Dxy in confined slit pores of HS governed by Newtonian 
dynamics®^. Details ofpg''® calculations are provided in 
the SM43. 

As is evident in Fig. 5, the average diffusivities D^y 
and Df'^{z) of the confined fluids approximately col¬ 
lapse onto the relevant bulk curves over many orders of 
magnitude in pg'"® for systems governed by either New¬ 
tonian and Brownian dynamics. Notably, even at high 
associated with supercooling, Dl''^{z) values only 
differ from the bulk by factors of 2-3 based on com¬ 
ponent Pg'"®; if one instead plots diffusivities against a 
less “microscopic” static property, e.g., component 



FIG. 5. (color online). Average component diffusivities 
in the xy-plane (squares) and z-direction (triangles) versus 
average component insertion probabilities pg''® for small and 
large particles (filled and unfilled symbols, respectively) cal¬ 
culated from MD simulations (main) and BD simulations 
(inset) of pore size H = 5 and average total packing frac¬ 
tions <)>^''® = 0.20, 0.30, 0.35, 0.40, 0.45, 0.48, 0.50, 0.51, and 
0.52. Average diffusivities for bulk mixtures shown as solid 
black lines with ±20% bounds shown as dashed lines. 


confined and bulk diffusivities differ by up to an order of 
magnitude. Overall, the data support the idea that, de¬ 
spite the difficulty of rationalizing position-dependent 
diffusivity behaviors based on bulk physics, the aver¬ 
age dynamics of inhomogeneous fluids are nonetheless 
strongly encoded with bulk correlations between mobil¬ 
ity and available space. 

In closing, by characterizing the particle dynamics of 
highly confined binary HS mixtures in both inhomo¬ 
geneous and isotropic dimensions, we find that diffu¬ 
sive mobility is not universally predicated upon packing 
structure according to bulk HS behaviors, as exempli¬ 
fied by the reversal from positive to negative correlations 
between local total packing fraction (j}{z) and single¬ 
particle diffusivity D^''^{z) coinciding with the onset 
of supercooling. In contrast, average diffusive mobil¬ 
ity is strongly encoded by the bulk physics, and can be 
approximately predicted via knowledge of the distribu¬ 
tion of available space. For the confined fluids studied 
here, results are insensitive to whether Newtonian or 
Brownian (i.e., overdamped Langevin) microscopic dy¬ 
namics govern particle trajectories, though it is an open 
question as to whether similar classes of behavior will 
emerge in real colloidal thin films treated within the 
FP formalism. More speculatively, the shapes of the 
D^''^{z) profiles (and their qualitative insensitivity to 
^avg) suggegt tJiat there may simply be a “universal” 
oscillatory signature of local diffusivity that emerges for 
non-continuum fluids proximal to confining potentials 
regardless of microscopic dynamics or the specific na¬ 
ture of any emergent structural inhomogeneity. We are 
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presently investigating this possibility. 
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APPENDIX A: SIMULATION PROTOCOLS 

To generate particle trajectories governed by either 
molecular dynamics (MD) or Brownian dynamics (BD), 
we simulate systems of Ngm + Mg = 2400 particles using 
GROMACS 4.5.5"^^. MD trajectories are generated by 
integrating the Newtonian equations of motion with a 
time step of 0.001 while fixing temperature with a Nose- 
Hoover thermostat. BD trajectories are generated via 
the overdamped Langevin equation (ignoring hydrody¬ 
namic interactions), where the position of particle i 
is propogated with a time-step of 0.01 according to®°’^^: 
r,{t + At) = Yi{t) + D^.AfF^{v^{t)) + Here, 
is the infinite dilution diffusivity, Fi{t) is the net force 
due to interparticle and wall interactions, and ^i(t) is 
the stochastic contribution. We set = 0.001 and 
= tTsm/o’ig, and in each direction, ^i{t) = 
r^{t) where r^{t) is a Gaussian noise with 

(r^(t)) = 0 and variance = 1. 

To generate bulk and confined packings at high 
we initialize systems at < 0.30 and compress 

them to the desired packing fractions via the method 
of Lubachevsky and Stillinger®^, in which particle di¬ 
ameters are grown linearly with time according to the 
dimensionless growth rate F. We execute compressions 
via MD simulations with effective F < 1x10“®, which al¬ 
lows us to avoid generating partially jammed (i.e., non- 
equilibrated) structures for all presented Further 

equilibration and production runs (MD and BD) are 
then initialized with the final structures. 


APPENDIX B: DENSITY PROFILES ACROSS 
CONFINED PORES 

In the main text, we largely plot local total pack¬ 
ing fractions </)(z) to characterize fluid structure. How¬ 
ever, a few noteworthy aspects concerning fluid struc¬ 
ture are apparent by simultaneously considering indi¬ 
vidual component density profiles p(z), as shown here 
in Fig. 1. First, for select conditions (e.g., = 0.40), 

a given individual component p(z) profile can be out of 
phase with (j){z). This does not happen frequently, and 
given the characteristics of the mixture (e.g., composi¬ 
tion Xsm = 0.75, size ratio uig/crsm = 1-3) studied here. 




FIG. 6. (color online). Local total packing fractions 
<j){z) (line-symbols) and individual component densities p{z) 
(lines) from simulations of pore size 77 = 5 at various aver¬ 
age total packing fractions where results for small and 
large particles are plotted with darker and lighter curves, 
respectively. 


this behavior is only observed for small particle pro¬ 
files. This underlines the importance of characterizing 
the spatial distribution of both particle species in order 
to gain a complete picture of packing structure; oth¬ 
erwise, one may come to qualitatively incorrect conclu¬ 
sions about which regions of a fluid are densely packed. 
However, for the size ratio studied here, qualitative vari¬ 
ations in (f){z) can mostly be derived from knowledge 
only of the large particle p{z) profiles (as would likely 
be possible for all Xsm < 0.75). 

As discussed in relation to Figs. 2 and 5 in the main 
text, by comparing these p{z) profiles against local par¬ 
ticle diffusivities D^''^{z), it is evident that there is no 
(^av®-independent local correlation between component 
density and diffusive mobility. Likewise, the local avail¬ 
able space for particle insertion, as quantified by inser¬ 
tion probability Po{z) in these systems, does not univer¬ 
sally correlate with mobility because (see Section IV) 
Po)z) = p{z)/^, where ^ is the spatially-invariant com¬ 
ponent activity^^. Of course, this implies that other 
position-dependent measures of structural correlations 
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FIG. 7. (color online). Local total packing fractions (l){z) 
(top) and local diffusivities in the ^-direction (z) of small 
(middle) and large (bottom) particles calculated from MD 
simulations of pore size H — 5 and average total packing 
fractions = 0.20, 0.30, 0.35, 0.40, 0.45, 0.48, 0.50, 0.51, 
and 0.52. 


that are positively correlated with po{z), such as the lo¬ 
cal two-body excess entropy s^^^(z), would likewise ex¬ 
hibit no consistent correlation with diffusive mobility^®. 
Taken altogether, the results in the supplemental Fig. 
1 underline that idea that local structure does not con¬ 
trol, nor can be used to predict, local diffusive mobility 
in any straightforward way. 


APPENDIX C: COMPONENT LOCAL DIFFUSIVITIES 

In supplemental Fig. 2, we show local particle diffu¬ 
sivities for both small and large particles gov¬ 

erned by MD at the various conditions studied. As 
noted in the main text, it is evident that the 
profiles for the two types of particles have qualitatively 
similar shapes and exhibit the same insenstivity to the 
structural rearrangment that occurs upon the onset of 
supercooling at high Note that there is slight drift 
(i.e., overestimaton of Dl''^[z)) near the edges of the 
large-particle profiles at higher due to the relatively 
meager amount of particle trajectory data obtained for 
this species (xig = 0.25). Df'^{z) prohles for small and 
large particles governed by overdamped Langevin dy¬ 
namics are also qualitatively similar. 


APPENDIX D: CALCULATING AVERAGE INSERTION 
PROBABILITY 


To calculate Pq'"® for each component, we note that if 
there is no external field at position z (i.e., = 0), 


the local insertion probability®^ for bulk or inhomoge¬ 
neous HS is a ratio"^’^ Po{z) = p{z)/^ of the local com¬ 
ponent density p{z) and the spatially-invariant compo¬ 
nent activity ^ = exp(/3p)/A®, where the latter is de¬ 
fined by the component chemical potential fi and the 
de Broglie wavelength A. Given that we have com¬ 
ponent p{z) profiles measured from the MD and BD 
simulations, all that is required to obtain po{z) pro¬ 
files are activities ^ for the bulk and confined mix¬ 
tures at the various values. These ^ values are 
obtained via grand canonical transition matrix Monte 
Carlo (GC-TMMC) simulations®^, with implementation 
details presented elsewhere®®. It is then straightforward 
to calculate Pq''® = po{z)dz, which for bulk 

mixtures is simply Pq''® = p/^. 


^R. Kjellander and S. Sarman, J. Chem. Soc., Faraday Trans. 
87, 1869 (1991). 

^J. N. Israelachvili, Intermolecular and Surface Forces (Aca¬ 
demic Press, New York, NY, USA, 2011). 

^J. Wu and Z. Li, Annu. Rev. Phys. Chem. 58, 85 (2007). 

^R. Roth, J. Phys.: Condens. Matter 22, 063102 (2010). 

®V. Krakoviack, J. Phys.: Condens. Matter 17, S3565 (2005). 
®G. Biroli, J.-P. Bouchaud, K. Miyazaki, and D. R. Reichman, 
Phys. Rev. Lett. 97, 195701 (2006). 

^A. J. Archer, J. Phys.: Condens. Matter 18, 5617 (2006). 

®S. Lang, V. Bo^an, M. Oettel, D. Hajnal, T. Franosch, and 
R. Schilling, Phys. Rev. Lett. 105, 125701 (2010). 

®W. Olivares-Rivas, P. J. Colmenares, and F. Lopez, J. Chem. 
Phys. 139, 074103 (2013). 

^'^S. Lang, R. Schilling, V. Krakoviack, and T. Franosch, Phys. 
Rev. E 86, 021502 (2012). 

^^S. Lang and T. Franosch, Phys. Rev. E 89, 062122 (2014). 
Hummer, New J. Phys. 7, 34 (2005). 

Mittal, T. M. Truskett, J. R. Errington, and G. Hummer, 
Phys. Rev. Lett. 100, 145901 (2008). 

^^C. R. Nugent, K. V. Edmond, H. N. Patel, and E. R. Weeks, 
Phys. Rev. Lett. 99, 025702 (2007). 

^®M. Hlnczewski, Y. von Hansen, J. Dzubiella, and R. R. Netz, 
J. Chem. Phys. 132, 245103 (2010). 

Carmer, F. van Swol, and T. M. Truskett, J. Chem. Phys. 
141, 046101 (2014). 

I'^Y. Rosenfeld, Phys. Rev. A 15, 2545 (1977). 

1®Y. Rosenfeld, J. Phys.: Condens. Matter 11, 5415 (1999). 

J. Hoyt, M. Asta, and B. Sadigh, Phys. Rev. Lett. 85, 594 

( 2000 ). 

^®R. Sharma, S. N. Chakraborty, and C. Chakravarty, J. Chem. 
Phys. 125, 204501 (2006). 

Mittal, J. R. Errington, and T. M. Truskett, J. Phys. Chem. 
B 110, 18147 (2006). 

Mittal, J. R. Errington, and T. M. Truskett, J. Chem. Phys. 
125, 076102 (2006). 

^^N. Gnan, T. B. Schrpder, U. R. Pedersen, N. P. Bailey, and 
J. C. Dyre, J. Chem. Phys. 131, 234504 (2009). 

24W. P. Krekelberg, M. J. Pond, G. Goel, V. K. Shen, J. R. 

Errington, and T. M. Truskett, Phys. Rev. E 80, 061205 (2009). 
^®M. J. Pond, J. R. Errington, and T. M. Truskett, J. Chem. 
Phys. 134, 081101 (2011). 

Nayar and C. Chakravarty, Phys. Chem. Chem. Phys. 15, 
14162 (2013). 

2'^J. C. Dyre, J. Phys. Chem. B 118, 10007 (2014). 

J. Mittal, J. R. Errington, and T. M. Truskett, Phys. Rev. Lett. 
96, 177804 (2006). 

J. Mittal, J. R. Errington, and T. M. Truskett, J. Phys. Chem. 
B 111, 10054 (2007). 

Mittal, V. K. Shen, J. R. Errington, and T. M. Truskett, J. 
Chem. Phys. 127, 154513 (2007). 

®^G. Goel, W. P. Krekelberg, J. R. Errington, and T. M. Truskett, 
Phys. Rev. Lett. 100, 106001 (2008). 


























7 


32g. Goel, W. P. Krekelberg, M. J. Pond, J. Mittal, V. K. Shen, 
J. R. Errington, and T. M. Truskett, J. Stat. Mech.: Theory 
Exp. 2009, P04006 (2009). 

Chopra, T. M. Truskett, and J. R. Errington, Phys. Rev. E 
82, 041201 (2010). 

^■*3. J. Borah, P. K. Maiti, C. Chakravarty, and S. Yashonath, 
J. Chem. Phys. 136, 174510 (2012). 

®®X. Ma, W. Chen, Z. Wang, Y. Peng, Y. Han, and P. Tong, 
Phys. Rev. Lett. 110, 078302 (2013). 

®®Y. Liu, J. Fu, and J. Wu, Langmuir 29, 12997 (2013). 

®^T. S. Ingebrigtsen, J. R. Errington, T. M. Truskett, and J. C. 

Dyre, Phys. Rev. Lett. Ill, 235901 (2013). 

®®K. V. Edmond, C. R. Nugent, and E. R. Weeks, Phys. Rev. E 
85, 041401 (2012). 

A. Bollinger, A. Jain, and T. M. Truskett, Langmuir 30, 
8247 (2014). 

A. Bollinger, A. Jain, and T. M. Truskett, J. Phys. Chem. 
B (2014), DOI: 10.1021/jp508887r. 

Chandler, J. D. Weeks, and H. C. Andersen, 220, 787 
(1983). 

^^B. Hess, C. Kutzner, D. van der Spoel, and E. Lindahl, J. 
Chem. Theory Comput. 4, 435 (2008). 


■^^See supplemental material at [URL will be inserted by AIP] for 
detailed simulation protocols and auxiliary data. 

Liu, E. Harder, and B. J. Berne, J. Phys. Chem. B 108, 
6595 (2004). 

■**0. H. Weiss, “First passage time problems in chemical physics,” 
in Adv. Chem. Phys. (John Wiley & Sons, Inc., 2007) pp. 1-18. 

Sedlmeier, Y. Hansen, L. Mengyu, D. Horinek, and R. Netz, 
J. Stat. Phys. 145, 240 (2011). 

Widom, J. Chem. Phys. 39 (1963). 

■**3. Sastry, T. M. Truskett, P. G. Debenedetti, S. Torquato, and 
F. H. Stillinger, Mol. Phys. 95, 289 (1998). 

■^^E. R. Weeks and D. A. Weitz, Phys. Rev. Lett. 89, 095704 

( 2002 ). 

L. Ermak and J. A. McCammon, J. Ghem. Phys. 69, 1352 
(1978). 

®^M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids 
(Clarendon Press, New York, NY, USA, 1989). 

®^B. D. Lubachevsky and F. H. Stillinger, J. Stat. Phys. 60, 561 
(1990). 

J.-P. Hansen and 1. R. McDonald, Theory of Simple Liquids, 
3rd ed. (Academic Press, New York, NY, USA, 2006). 

. K. Shen and J. R. Errington, J. Chem. Phys. 122, 064508 
(2005) 



